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VALIDITY OF THE EXPECTED EULER CHARACTERISTIC 

HEURISTIC 

By Jonathan Taylor, Akimichi Takemura and Robert J. Abler 

Stanford University, University of Tokyo and Technion 

We study the accuracy of the expected Euler characteristic ap- 
proximation to the distribution of the maximum of a smooth, cen- 
tered, unit variance Gaussian process /. Using a point process repre- 
sentation of the error, vaUd for arbitrary smooth processes, we show 
that the error is in general exponentially smaller than any of the 
terms in the approximation. We also give a lower bound on this ex- 
ponential rate of decay in terms of the maximal variance of a family 
of Gaussian processes derived from the original process /. 



1. Introduction. In this paper, we study the expected Euler character- 
istic approximation to 

(1) p(snpfix)>u] 

where / is the restriction to M of /, a process on a manifold M, and 
M is an embedded piecewise submanifold of M. 

When the process / is Gaussian with zero mean and has unit variance, 
the expected Euler characteristic approximation is given by 

pf sup/(x)>n) =E(x(Mnrn^^,+oo))) 

(2) 

dimAf 



^ /:,(M)(27r)"(^+i)/2 r Hj{r)e-'"/^dr, 

where Hj is the jth Hermite polynomial, / ""^[^j +oo)) is the Eu- 

ler characteristic of the excursion of / above the level u and the Cj{M) 
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are the intrinsic volumes, or Lipschitz-Killing curvatures of the parameter 
space M [20, 22], measured with respect to a Riemannian metric induced 
by /, which is discussed below in Section 2. In [20] only the special case 
of finite Karhunen-Loeve processes (see below) was treated and in [22] the 
case of manifolds with smooth boundaries was treated. The result, proved 
in the generality described above, is to appear in [4]. With a slight abuse of 
notation, we have used P to denote our approximation to (1); that is, we 
are using P in the statistical sense of a "point estimator" and not as some 
alternative probability measure. ^ 

As noted above, the case when M is C°° and / is a centered, unit vari- 
ance, finite Karhunen-Loeve expansion process was studied in [20]. These 
assumptions imply that there exist a map (/9:M — > 5'(M") where ^(M") is 
the unit sphere in and a random vector i{uj) ~ N{0,lnxn) such that 

n 

f{x,uj) = ((/?(x),^(u;))Kn = ^^j{u:)Lpj{x). 

In this setting, without loss of generality, we can assume that M is an 
embedded submanifold of /^(M") and ip is just the inclusion map. Using 
the volume of tubes approach [18, 23], it was shown in [20] that if M is a 
piecewise C°° submanifold of S'(M") the error in the above approximation 
is bounded by 



(3) 



P( sup f{x) > u 

c 

< 



P ( sup /(x) > u 



r(n/2)2("-2)/2 J^/^^,e^^M) 

= Cx P{xl > u^/ cos^ OciM)) 

where 9c{M) is a geometric quantity known as the critical radius of M 
[11, 12, 18, 21]. For M C 5(M"), the critical radius is roughly defined by the 



following property: for every < ^? < 9c{M) and arbitrary z € 

ds{Rn)(^,M)<0 = 



argmin(i5(iRn)(2;,x) is unique 



where the metric on 



is the geodesic metric 

ds{w^){x,z) = cos~'^{{x,y)). 

In another setting, when / is "almost" isotropic on M'^, then, with some 
additional assumptions on M (cf. Theorem 4.5.2 in [3]) Piterbarg [17] showed 
using the "double-sum" method that the error in using the expected Euler 
characteristic approximation is bounded by 

(4) P ( sup /(x) >u\-Y>\ sup fix) > < C7e~°"'/2 



P ( sup /(x) > u 
~xeM 
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for some a > 1, though no expression for a is given. In the one-dimensional 
case, more was previously known [5, 13, 16] (see below). 

Note that both bounds show that the error in approximating (1) is ex- 
ponentially smaller than all terms in the expected Euler characteristic ap- 
proximation. While undeniably useful, these two situations do not cover all 
possibilities. Referring to (4), certainly not every smooth Gaussian process 
of interest is isotropic, nor are the conditions required of M easily inter- 
pretable (cf. Definition 4.5.1 in [3]). Referring to (3), while every Gaussian 
process does admit an infinite orthogonal expansion, see [2] 

oo 

(5) /(x,a;) = ^Cj(u;)99j(x) 

i=i 

through its reproducing kernel Hilbert space (RKHS), it is clear that sub- 
stituting 71 = oo into (3) is meaningless. In fact, the situation is even worse 
in that the two cases do not even overlap: an isotropic field restricted to a 
bounded domain T C M^' cannot have a finite Karhunen-Loeve expansion 
[19]! 

This brings us to the main result of this work, Theorem 4.3, which, when 
/ is a constant variance Gaussian process as in the works cited above, pro- 
vides bounds for the error in using the expected Euler characteristic approx- 
imation to (1). Specifically, when / has unit variance, we show that 

lim inf —u~'^ log 

(6) 

Above, the "critical variance" cr'^{f) depends on the variance of an auxiliary 
family of Gaussian processes {f^)xeM, defined in (27) below. 

An alternative approximation to (1) is to use the expected number of 
(extended outward) local maxima [10]. The term "approximation" is used 
in a somewhat loose sense, as, to the authors' knowledge, there are no gen- 
erally applicable known closed form expressions for the number of extended 
outward local maxima of a smooth process. The only results known are 
ones which relate the asymptotic behavior of the expected number of local 
maxima of a Gaussian field on a manifold without boundary (which renders 
the qualifier "extended outward" unnecessary) to the expected Euler char- 
acteristic approximation, see [10]. Nevertheless, if one could compute the 
expected number of local maxima exactly, as one can the expected Euler 
characteristic in certain cases, one might expect to get a better approxima- 
tion to (1). Virtually identical arguments to those used in this paper show 
that when M is a manifold without boundary and / is Gaussian with con- 
stant variance, on an exponential scale the errors in the approximations are 
equivalent, though, in the interest of brevity, we do not pursue this here. 
When the manifold M has a boundary, the situation is more subtle and it 



P ( sup f{x) > u 



P ( sup f{x) > u 

.xGM 



1 1 

-2 + 2^' 
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may indeed be the case that the expected number of extended outward local 
maxima may be more accurate on an exponential scale. 

The critical variance is closely related to the critical radius appear- 

ing in (3). Specifically, when / is a centered, unit variance, finite Karhunen- 
Loeve expansion process, and M is a manifold without boundary, it is proven 
in Lemma 5.1 that 

al{f) = cot^e,{M) 

where Oc{M) is the critical radius of M, mentioned above and used in [21]. 

We note that, while (3) is an explicit bound for finite Karhunen-Loeve 
expansion Gaussian processes, it is not sharp, nor generally applicable. In 
particular it depends on the (generally unknown) dimension of the sphere 
into which M is embedded, that is, the dimension of the sphere in which 
M sits. In a companion paper [14], when / is a centered, unit variance, 
finite Karhunen-Loeve expansion Gaussian process, the asymptotic error 
as n ^ oo is evaluated using a Laplace approximation, rather than just 
the exponential behavior, which is the topic of this paper. In some one- 
dimensional stationary cases, the exact asymptotics of the error were found 
by Piterbarg in [16]. The results were generalized in [5], enlarging the class 
of processes covered by Piterbarg's result. Roughly speaking, the results of 
[5, 16] hold when the critical variance <Tc{f) the stationary process / is 
achieved locally. In this case, the critical variance is explicitly computable 
in terms of the spectral moments of / (see Lemma 5.3 below). 

Our main result. Theorem 4.3, is formally an application of Theorem 3.3 
to the case when / is Gaussian with constant variance. Theorem 3.3 gives 
a bound for the error of the expected Euler characteristic approximation 
for the restriction of an arbitrary suitably regular (cf. [1, 22]) process / 
on a manifold M to any embedded piecewise submanifold M C M. 
Theorem 3.3 is, to the authors' knowledge, the only available bound for the 
error in the expected Euler characteristic approximation for arbitrary, suit- 
ably regular, smooth random fields and should prove useful in studying the 
accuracy of the Euler characteristic approximation to non-Gaussian fields 
[8, 9, 24]. The analogy to (3) for nonconstant variance Gaussian fields using 
a variant of the volume of tubes approach is presented in [15]. 

Another noteworthy feature to our approach is that it is a direct approach 
to determining the error in using the Euler characteristic approximation. 
This should be contrasted with the bounds (3) and (4) which were both 
arrived at indirectly in the sense that the bounds were derived for the "vol- 
ume of tubes" approach and the "double-sum" approach and subsequently 
shown to hold for the expected Euler characteristic approximation. 

The proof of Theorem 3.3 depends on a point set representation for the 
global maximizers of /i = /i[Af , the restriction to M of a smooth deterministic 
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function /i : M — > M, above the level u. Of course, there is a trivial, and not 
very useful, point set representation of the set of maximizers of h above the 
level u: 



Lemma 2.2 gives an alternative point set representation of the maximizers 
using an auxiliary family of functions {h^)x£M- Once we have a point set 
representation of the maximizers of a smooth function, we apply a "meta- 
theorem" for the density of point processes arising from smooth processes 
of [1, 3, 4, 7, 22]. If the global maximizer of the process / is almost surely 
unique, the total mass of the density of the point process of maximizers above 
the level u is therefore just (1). The auxiliary family of processes {f^)xeM 
mentioned above is, in some sense, the stochastic analogue of {h^)x<=Ai in 
the deterministic setting. 

After defining the processes and describing their properties, we derive 
the following almost sure bound: 



where extended outward critical points are defined in a suitable fashion 
(cf. Section 2) and, for each x € M, the process is uncorrelated with 
f{x). Therefore, the points that contribute to the error are points where a 
Gaussian random variable f{x) and the supremum of a process f^{y), inde- 
pendent of f{x), are above the level u. The variance of the processes is 
what establishes the exponentially small relative error in (6). The only point 
where the Gaussian assumption is used is in bounding the expected number 
of points above, and the argument used here can be expected to extend to 
non-Gaussian processes as well. 

The organization of the paper follows. Sections 2 lays out the regularity 
conditions needed for Theorem 3.3, and reviews some notions of piecewise 
smooth manifolds. Theorem 3.3 is proved in Section 3, and Section 4 deals 
with the unit variance Gaussian case, where Theorem 4.3 is proven. We 
conclude in Section 5 with some examples; specifically we compute 
for stationary processes on M and isotropic processes on M'^ restricted to 
compact convex subsets. 



X G M : h{x) = maxh{y),h{x) > u 



^snp,^„ f(x)>u} - X{M n f \u,+00))\ 



<#<a;GM:xisan extended outward critical point of / 




6 



J. TAYLOR, A. TAKEMURA AND R. J. ABLER 



2. Suitably regular processes on piecewise manifolds. In this section 
we describe the class of processes to which Theorem 3.3 will apply. Before 
setting out our assumptions, we recall some basic facts about piecewise 
{I > 2) submanifolds of an ambient (j > manifold M. For a more 
detailed self-contained treatment of the material in this section, specifically 
as it relates to smooth processes, the reader is referred to [4]. 

A fc-dimensional piecewise manifold has a (not necessarily unique) 
decomposition into i-dimensional manifolds without boundary, <i <k; 
namely, 

k 

M=\J dMi. 

i=0 

Associated to every point x € M is its support cone in T^M , the tangent 
space to M at x 

S^M = {X^(^T^M:35>^,ceC^{{-5,5),M), 

c{t)eM VtG [0,(5),c(0) =x,c(0) = X^}. 

In words, the support cone of M at x is the set of all directions in which 
a smooth curve can leave x into M, but, in an infinitesimally small time 
period, still remain in M. 

When M is endowed with a Riemannian metric it is possible to define the 
dual cone of SxM . The following Riemannian metric will be essential to our 
ana^sis: an differentiable process induces a natural Riemannian metric 
on M given by (cf. [22]) 

g^iX^.Y,) = Coy{X J,yJ), 

forallX,,y,Gr,M. _ 
The dual cone of SxM, SxM* , in this case called the normal cone in M 
at X G M, is defined by 

NxM = SxM* = {Xx G TxM :gx{Xx,Yx) < 0, VF,. G SxM}. 

For X G dMk = M°, NxM = {TxM°)^, the orthogonal complement of dMk 
in TxM. The normal cone figures prominently in the approximation, and in 
the main result of this paper, as both the Euler characteristic and the global 
maximizer point processes are defined in terms of extended outward critical 
points, that is, critical points at which the gradient (viewed as a tangent 
vector in the ambient space) is in the normal cone of the set at the critical 
point. Roughly speaking, this means that the function is increasing along 
curves leaving the set along certain "normal directions." 

To avoid trivialities, we further assume our decomposition of M is such 
that for every <i <k and each x G dMi , the tangent space dMi is the 
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largest subspace contained in SxM . This condition rules out trivial decom- 
positions of a nice open set O C into {O \ F,F} for a finite point set or 
some smooth closed curve F cO. This condition is not strictly necessary for 
the Morse theorem of [20]. In its place, we could impose a condition on the 
Morse functions to not have any extended outward critical points on such 
F's. As this will almost surely be the case for the paths of suitably regular 
processes this point is somewhat moot. 

For our purposes here, piecewise smooth manifolds are required to have 
the additional property that they are locally approximated by SxM in the 
sense that, for every x G M, there exists a diffeomorphism 

(7) ^x:UxCM^MUx)cTxM 

of some neighborhood Ux of x ^ M, whose inverse, when restricted to SxM, 
maps any sufficiently small neighborhood of the origin to a neighborhood of 
X G M. This condition, for instance, rules out cusps in M as can occur when 
two manifolds intersect nontransversally. The importance of this condition, 
besides the fact that it is necessary in order to use the Morse theorem of 
[20], is that it implies that the set 

(8) {Xx G N{dMi) : Xx G {P^^qm^NxMY} 

is open in N{dMi) where {Pj^^Qjyj,NxM)° is the (relative) interior oi Pj^^qj^.^.NxM 

in TxdM^ and N{dMi) is the normal bundle of dMi in M. 

Having defined piecewise smooth manifolds, we now set out the assump- 
tions on the processes / defined on our piecewise smooth manifolds. 

Assumption 2.1. We assume that^ is the restriction of / to M, where 
/ is a square-integrable process on M, a g-dimensional manifold, and 
M is a compact, embedded piecewise /c-dimensional submanifold of M 
such that SxM is a convex cone for each x G M. We further assume that, 
for each i, the gradient of f\QMi read off in some nonrandom orthonormal 
frame field Ei = (Xi, . . . ,Xi) on dMi satisfies the conditions of Lemma 2.5 
of [22]. We denote this M*-valued process by V f\QMi,Ei - These conditions are 

satisfied if f\QMi is suitably regular in the sense of [22]. Finally, we assume 
that 

p{x,y) = CoY(f{x)J{y)) = l <^ x = y. 

At this point, it is probably worth describing the significance of the above 
assumptions. In point form, the significance of these assumptions are: 

1. The requirements that f = f\M foi^ ^ process on M and that the 
support cone SxM is convex for each x G M are necessary for the point 
process representation of the Euler characteristic and to apply the Morse 
theorem of [20]. 
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2. The conditions on the process f\dMi,Ei ensure that the expected number 
of critical points of f\QMi is finite and the density of point processes based 
on the critical points has an explicit integral representation in terms of a 
point process "meta-theorem" [3, 22]. 

3. The condition on p{x,y) ensures that the map x i-^ f{x) is an embedding 
of M into L^{n, T, P) which rules out "global" singularities in the process. 



2.1. A point set representation for the global maximizers of deterministic 
functions. Although our primary interest is with stochastic processes, we 
begin here with the deterministic case, stating and proving Lemma 2.2 which 
describes a point process representation of the global maximizers of a de- 
terministic function h. However, along each stratum dMi the point process 
depends on a function defined on 

dMi xM~ = dAk X M \ {(x, y) e dMi x M:x = y} 

which is singular near the "diagonal" {{x,y) G dMi x M ix = y}. We resolve 
this singularity in Corollary 2.6, and arrive at the point process representa- 
tion for the global maximizers in Corollary 2.7. 

Lemma 2.2. Suppose h = h^^j, the restriction of h £ C'^{M). Fix x G 
dMi, and choose € C'^{M, (— oo, 1]) such that 

a'^{y) = l =^ h{x)=h{y), a''{x) = l. 

Then, x is a maximizer of h above the level u if, and only if, the following 
three conditions hold: 

(i) h{x) > u. 

(ii) Vh[x) G NxM . That is, x is an extended outward critical point of h. 

(iii) h{x) > supj^gA/\{2,} h^{y), where, for all y ^M 

i/a^(y) = l. 

Further, ifV'^a^{x) is nondegenerate, and x is a critical point of h\Qj^., 
then, for any curve c:{—5,5) dMi with c(0) = x,c(0) = Xx ^T^dMi, 

V=/i|8„,(i)(.Y„A'J-V'af(,„(i)(X„XJft(j,) 
= -V^a|.,(.)(k.X,) • 

Remark. The condition that a^{x) = 1 ensures that for each x G M, x 
is a critical point of . 
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Proof of Lemma 2.2. The condition h{x) > n is self-evident. Suppose, 
then, that x E (9Mj,0 < f < A;, is a maximizer of h. Then V/i(x) G N^M; 
otherwise there exists a direction Xx € SxM such that gx{Xx,Vh{x)) > 
and X cannot be a maximizer. 

Because x is a maximizer, for all y such that a^{y) < 1 it follows that 

h{y) -a''{y)h{x) 



1 - a^(y) 



< h{x). 



On the other hand, if a^{y) = 1, then, by choice of a^, h{y) = h{x) which 
proves that 

h{x) > sup h^{y). 

yeAI\{x} 

To prove the reverse implication, assume that x is an extended outward 
critical point of h^gj^j. and 

h{x) > sup h''{y). 
y€M\{x} 

Now suppose that x is not a maximizer of h; then there exists y £ M \ {x} 
such that 

h{x) < h{y). 

In particular, for such a y, our choice of implies that a.^{y) < 1. It follows 
that 

h{y)-a''{y)h{x) 

which is a contradiction. 

The limit (9) follows from two applications of I'Hopital's rule. Specifically, 
we note that x is a critical point of h\QMi by assumption and the properties of 

imply that it must also be a nondegenerate critical point of . Therefore, 

h{c{t))-a'{c{t))h{x) _ {d/dt){h{c{t))-a''{c{t))h{x)) 
t™ 1 -a^(c(t)) ~ d/dt{l-a''{c{t))) 

_ {d^ /dt^){h{c{t)) - a''{c{t))h{x)) 
"t™ {d'^/dt'^){l-a'^{c{t))) ■ 

Since x is also a nondegenerate critical point of 1 — a^(-) the conclusion 
now follows from the fact that, for any /3 G C^((9Mj) and any curve 
c : (-(5, 8) dMi with x = c(0) a critical point of /3 G C^{dMi) and c(0) = 

lim^/?(c(t)) = v2/3(x)(X,,X,). □ 
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Condition (iii) above will be crucial to our later results as it is the condi- 
tion which determines whether a given critical point is indeed a maximizer 
of h. The condition is not quite "ready to use," as we will need to consider 
an analogy to the quantity 

(10) W{x)= sup 

y&M\{x} 

as a function of x. However, it is not hard to show, again by two applications 
of I'Hopital's rule, that 

(11) Vh{x)(^N^M <=^ sup h'^iy) <oo, 

yeM\{x} 

(12) Vh{x)£N^M =^ inf /i^(y) = -oo. 

yeM\{x} 

In other words, the above quantity is only finite at extended outward critical 
points of h and the process is singular even if '\7h{x) € NxM . 

Although neither Lemma 2.2 nor (11) is exactly what we will need later, 
they contain a somewhat simplified version of our later arguments. 

2.2. Continuity of W{x): the ^^blow-up^^ of h. As previously noted, the 
function 

(13) (:c,y)^/i^(y) 

is singular when y is near x even though it is on dMi x M~ . In particular 

lim h"^ (y) 

y^x 

is undefined. Ultimately, for the point process representation, we are inter- 
ested in the continuity of W{x) on dMi, and this singularity can make the 
arguments a little delicate. ^ 

In this section our main result shows that W{x) is continuous. Our strat- 
egy is to "blow up" M around a neighborhood of x G 5Mj, and to relate 
h^{y), for X in this neighborhood, to a continuous function on the "blow- 
up." We use the term "blow-up" as we relate the function h^{y) to a "desin- 
gularized" version of h on the blow-up of M near x. 

By assumption, the parameter set M is locally approximated by its sup- 
port cone and we can, without loss of generality, describe the blow-up of M 
under the assumption that M C M = IR"^ and for some neighborhood Ux 

f/^ n M = {x} e SxM n u^. 

We can assume this because we have assumed in (7) that every x € M 
has such a neighborhood, and to establish continuity of VF(x) only local 
properties (in x) of h^{y) are needed. 
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For X G dMi, the support cone SxM contains the z-dimensional tangent 
space Tx dMi. By our assumptions on the decomposition of M, this subspace 
is maximal in SxM. Therefore, we can decompose SxM as 

Tx dMi X Kx 

where Kx is a convex cone that contains no subspace. 

For ease of exposition, for the moment we move to a simple Euclidean 
setting and replace SxM by 

K = Lx KcR'i 

with L a subspace, and K a convex cone containing no subspaces. We will 
describe the "blow-up" construction first in this scenario, and then move 
back to piecewise smooth spaces. In what follows, h = will be the re- 
striction of a generic function on to Wi. We first define the "blow-up" 
B[K) of K along L as the disjoint union of the spaces 

^ = Lx(i^\{0}), 
dX = Lx {Kr\S{W)). 
Above, X should be thought of as the image of 

L X K\{{x,y) ^ L X K -.x = y}, 

under the map 
(14) 

and the second space dX as the "boundary" of X . The boundary is attached 
as follows: a sequence of points {xn,yn) G X converges to (xo,yo) £ dX if, 
and only if, 

Vn 

\\yn\\ 

This notion of convergence corresponds to a sequence "i^~^{xn, y-n) = {xn,Xn + 
Un) converging to the diagonal 

{{x,y)eLx K:x = y} 

along a well-defined direction. 

Remark. Identifying the tangent bundle r(M^) with W x M'?, we can 
think of B[K) as a subset of TiMfl), in which the y's above are replaced 
with Xx S SxM C TxM3. When convenient, we will consider B{K) as either 
a subset of M^? or a subset of T{W). 
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In Lemma 2.2, {a^)x^M was a family of functions which did not necessarily 
arise as the partial map of a function a: M x M ^ (—00, 1]. In some cases, 
particularly in the stochastic setting below, this is a natural assumption to 
make. 

Lemma 2.3. Suppose a € C'^{W^ x M'') is such that the partial map 

o^^{y) = a{x,y) 

satisfies the conditions of Lemma 2.2 at every x £ L and such that the 
Hessian of the partial map is nondegenerate at every x L. Then, any 

h G C^(M'^) maps to a continuous function on B{K) as follows: 

^aKr \ h{x + y)-a{x,x + y)h{x)- {Vhx,y)ui 



h-'''{x,y) 



1 - a{x,x + y) 
and for (x, y) G dX 

hm h-^^{t,s) = h{x)+ 



{t,s)^{x,y) ' -V^a{x){y,y)' 

Proof. Two applications of rHopital's rule. □ 
The term 

(V/ia;,y)iR5 



1 - a{x,x + y) 

above "resolves" the singularity along the diagonal in some sense. In effect, 
it forces every x G L to be a critical point of the map 

h{x + y) - a{x, x + y)- h{x) - lyh{x),y)M.i 



m3 3ys 



1 — a{x^ x + y) 



Our motivation for introducing /i"'^^ is to describe the singularities in the 
function h^{y) at critical points x of h\L (recall from above that L takes 

the place of T^dMi in our general cone K). We are therefore interested in 
points where a; is a critical point of h^i = h^i . Note that if x is a critical 
point of then, for all y G M'^ 

h{x + y)-a{x,x + y)-h{x)-{Vh{x),y)ui 
h {x + y) = 



(15) + 



1 — a{x, x + y) 

{P^Vh{x),y)^, 
1 — a{x,x + y) 



h'^^^{x,y) + ^^tl]}i^^ 
1 - a[x,x + y) 
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where represents orthogonal projection onto L"*", the orthogonal com- 
plement of L in M3. The expression (15) indicates that at critical points x 

of h\i^, h^iy) is the sum of a well-behaved term, /i"'^ and a singular term. 
The above relation holds for all critical points of However, for our later 
arguments we would like to have this relation hold for all x G L to avoid 
having to condition on a set of measure zero. We therefore redefine as 
follows: 

, f h{y) - Q!'{y)h{x) - {PlV h{x) , y)Ki ;f^x/, n/. 
(16) h^{y) ^ l-a-{y) ' " ^ ^' 

[h{y), ifa-(y) = l. 

With this redefinition of h, relation (15) holds for all a; € L and, for each 
critical point x of the two definitions of coincide. For the remainder 
of this section, we will use the definition (16). 

Lemma 2.4. IfPi;S/h{x) e K*, the dual of K C , then for any hounded 
neighborhood Ox 3 x 

sup h^{y) < +00. 
yeOAi^} 

Proof. If P^Vh{x) e K*, then (15) implies that for ah y 

/i^(y) </i"'^(x,x + y) 

and h°''^ is continuous in y, and therefore bounded on bounded sets. If, on 
the other hand, Pi^Vh{x) ^ K*, then there exists a unit vector v ^ K such 
that {Pi^Vh{x),v)]Rq > and x + tv G Ox for sufficiently small t. Relation 
(15) implies that 

\imh^(x + tv) = +00. 

This follows from the fact that the numerator in the expression on the right- 
hand side of (15) is strictly positive of order 0{t) for t small, while the 
denominator is of order O(t^). □ 

Lemma 2.5. If P^\/h{x) S {K*)° , then for any bounded neighborhood 
O of the origin in M.'^ 

Woix)= sup /i (y)= sup {h"''^{x,v)+ ^ — — ] 

y£Kn{{x}®o)\{x} v£Kn{O\{0})^ ^ a[x,x-\-v) J 

is continuous at x. 
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Proof. We first note that by (15) the two suprema above are equal, 
and it suffices to consider the supremum on the right. 
Consider a convergent sequence 

{x,Vnix))^>Q {x,V*{x)) 

in B{K) along which the supremum Wo{x) is approached. Then either 
||^n(a^)||K>j >0 for all n sufficiently large, or ||t'n(3;)llR9 ^0. In the ffist case 
it is immediate that (x,w*(x)) is in X and, as we will show in a moment, in 
the second case {x,v*{x)) G dX r\L x S{L), where S{L) is the unit sphere 
in L. In other words, the limiting direction v*[x) is in S{L). Further, if 
||^^n(ic)||K>j — > 0, then it is easy to see that the sequence {x,PLVn{x)) also 
achieves the supremum Wo{x). 

To see why v*{x) must be in S{L) in the second case, suppose that 
{x,Vn{x)) converges to {x,v*{x)) with v*{x) £K\L. As P^Vh{x) G {K*)° 
it follows that 

{ptvh{x)y{x))^,<Q. 

Therefore, applying the same argument as in (15), 

{PtVh{x),Vn{x))^, 

hm -— — = — cx), 

n^oo \ — a{X,X + Vn[x)) 

which contradicts the assumption that 

lim h^{x + Vn{x)) = Wo{x). 

n — >oo 

The above implies^that for any x S L, and any convergent sequence 
{x,Vn{x)) achieving Wo{x)^ we can either assume that ||t)n(a;)|| is bounded 
uniformly below, or that Vn{x) € L for all n. Continuity of Wo{x) now fol- 
lows, as, for such sequences and e > sufficiently small 

sup \ff^^{x,Vn{x)) - B(h,a,K){y,Vn{x))\ < e, 

y€B{x,5{e)) 



sup 

y&B(x,S(s)) 



{P^Vh{x),Vnix))u'^ {P^Vh{y),Vn{x)h, 



1 - a{x,x + Vnix)) 1 - a{y,y + Vn{x)) 



< e. 



□ 



Returning to piecewise smooth spaces, the above arguments will need 
slight modifications. Specifically, the map ^, defined in (14), has no natural 
replacement candidate for a piecewise smooth space. However, as noted in 
a remark above, we can think of B{K) as a subset of the tangent bundle 
r(M'^) in which case ^ -.W ^ T{W). In the piecewise smooth setting, 
we must therefore replace the map ^ with a map H -.M x M ^ T{M ) such 
that for each x G M, 

(17) H{{x} xM) ctM. 
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One of the key properties of ^ was that the sequence converges 
to a point in as n — > as long as ^ x and the unit vector 

Vn 



WHn — XnWm 

converges in S(M''). We can replace this property of ^ by asking the following 
of H. For any curve c:(-<5,<5) M x M with c(0) = (ci(0), C2(0)) = 
(x,x),C2(0) =Xx, we require that 

(18) ^■^g(x,C2(t))-ff(x,x)^ 
^ ^ t^o t 

Given such an H, we can, as in (16), redefine for a function on a 
piecewise smooth space M. For x G (9Mj, we redefine as follows: 

/t(;/) -a^(y)/i(x) -g(V/t|aM^(x),ij'(x,y) -a^(y)g(x,x)) 

fig) h''(v) = l ^-a^'iy) 
^ ^ ^ if«-(y)y^i, 

h{y), if a^(2/) = 1. 

With redefined, by working in suitably chosen charts, it is not difficult 
to prove the following. 

Corollary 2.6. Suppose that a e C^(M x M) is such that the partial 
map 

a^(y) = a{x,y) 

satisfies the conditions of Lemma 2.2 at every x £ M and such that the 
Hessian of the partial map is nondegenerate at every x G M. Further, 
suppose H : M X M ^ T{M) satisfies (17), (18) and let h^ be defined as in 

(19) . Then, the first conclusion of Lemma 2.2 holds if we replace (i) with: 
(i') h{x) -g{H{x,x),Vh\gM, {x)) > u. 

The second conclusion of Lemma 2.2 reads: for every x € dMi and any 
curve c: (-5,5) ^ dMi with c{0) = x,c{0) = X^, 

^^^^^^^^^^ _VHh{-) -giH-i.),Vh^9MAxmX,,X,) 



(20) -V2aj^g^,^^(x)(X^,X^) 

-{h{x)-g{H^x),Vh\9MA^))) 

where H^{-) is the partial map H^{y) = H{x,y). 
Furthermore, the function 

W{x) = sup h^{y) 
yeM\{x} 
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is, for each <i <k, continuous on the set 

{x G dMi : P^^Q,,Vh{x) e {P^^QM^N.Mn, 

where Pji^gMi represents projection onto the orthogonal complement ofT^dMi 
in T^M. 

2.3. A point process for the global maximizers of stochastic processes. For 
the remainder of this work, we choose a fixed piecewise smooth space M and 
a process f on M D M satisfying Assumption 2.1. 

In this section, we describe process analogies of /i^, a and H in the case 
for which h is replaced with the smooth process /. Specifically, we take 

a{x,y) =p{x,y), 

(21) ^ 1 

H{x,y) = F{x,y) ^ ^ Cov(/(y), 

for some orthonormal frame field {Xi^x, ■ ■ ■ iXq,x) on TxM. 

The following is the stochastic analogy of Lemma 2.6, that is, the point 
process representation of the maximizers of /. 

Corollary 2.7. Under Assumption 2.1, almost surely, the maximizers 
of f are isolated and the maximizers of f are the points x £ 9Mj, <i<k, 
such that: 

(i) V/|aM,(x)=0; 

(ii) f{x)-g{F-{x),Vf\QM,{x))>u; 

(iii) PlQj,jyf{x)e{Pl9j,,NxMY; 

(iv) f{x)>svvpyf.M\{x}F{y), where 

( f{y) - p{x,y)fix) -g{F-{y) - pix,y)F-{x),S/f^9M,ix)) 

(22) ny)^\ ,„ , 1-P(x,2/) 

ifp{x,y)^l, 

^f{y), ifp{x,y) = i- 

Note. If the joint density of \/f{x), read off in some orthonormal basis 
of TxM, is bounded by some constant K uniformly in x € M, then, almost 
surely, there will be no critical points of f\dMi such that P-r^QMi^ fi^) € 
dNxM C TxdM^-. Therefore, almost surely, all global maximizers will be 
such that P'^^gj^.'V f{x) is in the relative interior of NxM in TxdM^. As 
for the proof of Lemma 2.10 below, the proof of this is reasonably standard, 
and follows along the lines of similar results in, for example, Chapters 11 
and 12 of [4]. We therefore omit the details. 
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Proof of Corollary 2.7. The only part of the argument in Lemma 2.2 
that needs to be modified is what happens when a^{y) = p{x,y) = 1. In 
the deterministic case, we assumed that a^{y) = 1 implied h{x) = h{y). 
In the random case, we know that a^(y) = 1 implies f{x) — E(/(x)) = 
(/(y) - E(/(y)))cr(2;)/cr(?/) almost surely, where 

Almost surely, then, it is stih true that if x is a maximizer of /, then /(x) > 
f^{y) for all y such that p^{y) = 1, since otherwise x cannot be a maximizer. 
The reverse implication follows similarly. □ 

With / defined as above, rt is easy to see that Corollary 2.6 holds almost 
surely with replaced hy , a hy p and H hy F. 

2.4. Point process representation for the difference between the expected 
EC and the true probability. Our assumptions ahow us to use the Morse 
theorem of [20] to express the expected Euler characteristic of the excursions 
M n f~^[u, +oo) as integrals over M. The formula is not new, though we re- 
peat it here for use in deriving bounds on the error in the Euler characteristic 
approximation. 

What is new, and is crucial to the entire paper, is the exact expression in 
Proposition 2.9 for the supremum distribution (1). 

Proposition 2.8. Under Assumption 2.1, and with the notation there, 

P( sup f{x) > u 1 

= E(x(Mnrnn,+(^))) 
(23) ^ 

= ^ / E(det(-V2/|OA/„i?,(^))lAfHV/|9M„i?,(2:) = 0) 

where TLi is an i-dimensional Hausdorff measure induced by g, W/igj,/. ^ W 
is the density of V f^QMi^E^^) o.'^d 

A^^ = {fix) >u,Vf{x) e N.,M}. 

Suitable regularity of the process V f\QMi,Ei{^) implies that the maximiz- 
ers^f / are almost surely isolated, though it does not guarantee uniqueness. 
If W{x) were continuous when restricted to (9Mj, Assumption 2.1 would al- 
low us to apply the general point process Lemmas 2.4 and 2.5 of [22] to the 
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point process representation of the maximizers in Corollary 2.7. The almost 
sure analogy of Corollary 2.6 shows that W{x) is not continuous, but it is 
continuous on the open set 

{x e dMi ■ PisMyfix) G {N^MY}. 

Further, we are only interested in its behavior on this set. Straightforward 
modifications of the above cited lemmas, which we omit, lead to the following 
representation for the supremum distribution. 

Proposition 2.9. Suppose that, almost surely, f has a unique maxi- 
mum and that Assumption 2.1 holds. Furthermore suppose that, for every 

xeM, 

P{W{x) = f{x)) = F{W{x) = u) = 0. 

Then, 

P sup f{x) > u I 
V x&M J 

(24) =Y,I E(det(-VV|9M.,i^,(^))lAft'HV/|aM.,i^,(x) = 0) 

><'^v/|s„,,B^(x)(0)ciWi(2;), 

where 

= {f{x) >uy W{x), V/(x) G N^M}. 

The discrepancy between the expected Euler characteristic approximation 
and the true supremum distribution is 

J^\Sf^M{u) = P f sup fix) > n) - P f sup f{x) > u) 
KxeM J \xeM J 

(25) = E / E(det(-VV|9Af.,i?,(a:))lAf««|V/|9M„z,,(x) = 0) 

X'/'v/|a„^,s^(x)(0)d7^i(x) 

where 

A™« = {u< f{x) < W(x),Vj(x) G N^M}. 

Before concluding this section, we provide a lemma giving sufficient con- 
ditions for the uniqueness of the global maximum of /. As mentioned above, 
the proof is reasonably standard fare and so omitted. Detailed arguments 
for very similar cases can be found in Chapters 11 and 12 of [4]. In the 
nonmanifold setting, these arguments are classical. (E.g., see Theorem 3.2.1 
of [1] or [6] and references therein.) 
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Lemma 2.10. Suppose that for allpairs 1 < i,j < k and allpairs {{x,y) :x G 
dMi , y € dMj } the random vector 

V{x,y) = (fix) - fiy),Vfi9M,,E.X^),^f\dM„E,iy)^ 
has a density, bounded by some constant K{x,y). Then, 

P({3 (x, y):xedMi,ye dMj,V{x, y) = € M^+-'+i}) = 0. 

3. Bounding the error. Expression (25) is an explicit formula for the 
error in the expected Euler characteristic approximation. A similar explicit 
formula can be derived for the error of the approximation based on the 
expected number of local maxima above the level u, though^ we do not pursue 
this here. However, as described in Section 2, the process is singular near 
X, and actually has infinite variance near x, which means that standard 
tools such as the Borell-Tsirelson inequality cannot be used to bound its 
supremum distribution. _ 

To see that the process can have infinite variance, assume that / is the 
restriction of an isotropic field to [0, 1]^. Fix a point (x,0). We will compute 
the variance of along the curve 

c(t) = (x,0)-t-(l,0). 

Straightforward calculations show that 

F-(c(t)) = 0, F-(x)=0, 

so that 

Ja;^^^^^)^ ^ fjy)- p{x,y)fix) 
^-p{x,y) 

In this case, the variance of f^{y) is easily seen to be 

Var(r(.)) = i±4^ 
^- p{x,y) 

and 



limVar(r(y)) = lim- 



'i- + pix,y) 



+ 00. 



p{x,y) 

In general, if x € dMi, then, along a curve c : {—6, 0] — > M with c(0) = —X^ G 
S^M\T,dMi, 

limVar(/^(c(t))) = +oo. 

Although this is somewhat worrying, in (25) we only care about large 
positive values of and, further, we only care about the behavior of on 
the set 

(26) {ix,u;): P^^eM,Vf{x){Lo) G {P^^gM^N^Mf}. 
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We exploit these facts and introduce a process in this section which 
has, under some conditions, finite variance and dominates on the set (26). 
It is this process whose variance appears in the exponential bound for the 
behavior of Diff in the Gaussian case. 

Obviously, the process which we define below, does not dominate the 
absolute value of Indeed, if this were true, the process would have 
infinite variance as well. 

The process is defined as follows: 

(27) r(y) = \ l-P^(y) 

I f{y) - 9{F^{y) - PN.MF'-{y),Vf{x)), if p-{y) = 1, 

where Pn^m ■ T^M N^M represents orthogonal projection onto NxM and 
F:M xM ^ T{M ) is given by 

,A _ /^(a;,y) -p(a;,y)i^(x,x), ifp(j;,y)/l, 

where F is defined in (21). 

Lemma 3.1. On the set (26) of extended outward normal points, for 
every y G M , 

ny)>r{y). 

If X ^ dMk = M° , then equality holds above. 
Proof. First, we note that 

(1 - P'^iy)) ■ {f{y) - ny)) = diF'iy) - PN.MF^{y),P^^9Mym)- 

As NxM is a convex cone, it follows that for any S TxM 

Yx - Pn.mYx G NxM* 

where NxM* is the dual cone of NxM, which is just the convex hull of 5a;M. 
By duality, 

9iYx-PN^MYx,Vx)<0 

for every Vx G NxM. Consequently, on the set (26) 

g(.Yx - PN^MYx,PiaM,Vf{x)) < 
for every Yx S TxM. As F^{y) G TxM for each y, the first claim holds. 
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As for the second, if x G dMk, then N^M = T^ dM^ and Pn^mVx = for 
ah Vx G Tx dMk. Similarly, 

g{Vx,Pi8M,vKx)) = Q. 

Therefore, on this set 

g{F-{y) - PM^MF-{y),vKx)) =0. □ 



As far as the continuity (in x) of 

y&M\{x} 



(29) W{x)= sup f 



is concerned, it is not difficult to show that, almost surely. Corollary 2.6 
holds with replaced by that is, that W{x) is continuous on the set 

{x G dMi : PiaM^m e {P^QM^NxMY}. 

As we will see in the proof of Theorem 3.3, Lemma 3.1 provides the basic 
bounds for Diff j^jv/(^^)- The following corollary to Lemma 2.2 will also be of 
use to us. 

Corollary 3.2. // / has unit variance, then, for any unit speed 
curve c : {—6, 6) — > dM^ with c(0) = x, c(0) = Xx, 

V^fm,Ax){Xx,Xx) - V^p^{x){Xx,Xx)fix) 



-o^^(^(^»- -VV(x)(X.,X.) 

■|9Mfe( 



-Vy\eMAx){Xx,Xx) + f{x). 



Further, 

(30) sup |-VV|0A4(x)(X„X,) + /(x)|< sup \r{y)\. 

X^&S{T^dMk) S/eA/\{x} 

Proof. The proof is essentially just the second conclusion of Lemma 2.2, 
recast in the stochastic process framework. The only thing that needs to be 
verified is that 

V^p^{x){Xx^Xx) = -l, 
but this follows from the fact that 

V^p-{y){Xy,Yy) = CoY{V^f\QMMiXy,Yy)J{x)) 
and the fact that, double form 

Cov{V^f\QM,{x)[Xx,Yx),f{x)) = -gx{Xx,Yx), 

(cf. [22]). □ 

Using the results of Lemma 3.1 we have the following theorem. 
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Theorem 3.3. Suppose that f has a unique maximum, almost surely, 
and Assumption 2.1 holds. Further suppose that, for every x G M , 

P{W{x) = f{x)) = P{W{x) = n) = 0. 

Then, 

|Diff/,A/HI 

(31) <Y,f E{\det{-V^f\aMMx))\lsERn\Vf^eA.u,ESx) = 0) 
where 

B^^^ = {u<f{x)<W{x)}. 

4. Gaussian fields with constant variance. In this section, using Theo- 
rem 3.3, we derive an explicit bound for the exponential behavior of Diff /^^/(u) 
when / is a Gaussian field with constant variance, satisfying Assumption 2.1. 
The assumption of constant variance implies certain random variables are 
uncorrelated, hence independent in the Gaussian case. In particular, the 
assumption of constant variance implies that for x E 9Mi,0 <i<k, the en- 
tire process {f^{y))y^M\{x} is independent of f{x) as well as ^ f\aMi{x)- 
This allows us to remove the conditioning on ^ f\dMi{x) below. Once this 
conditioning is removed, the rest of the argument relies only on the Borell- 
Tsirelson inequality [2]. 

Our first observation is that, whether / has constant variance or not, 
for each x € M, the process f^{y) is uncorrelated with the random vector 
V f\dMi{x)- Hence, in the Gaussian case, f^{y) is independent of V/|gjvf^(x). 

Lemma 4.1. For every x € 9Mj,0 <i<k and every y ^M \ {x} 

Cov(r(2/),X,/) = 

for every € dMi . 

Proof. We first note that, if p^{y) / 1, then 
(l-p-(y))Cov(r(y),X,/) 

= Cov(/(y)-p^(y)/(x),X,/) 

- Cov(^(i?-(y) - P;v.mF-(2/), V/(x)),X,/) 
= Qo-^{f{y)-p^[y)f{x),X.,f)-g{F'\y)-P^^MF^{y),X^,). 
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If, on the other hand, p^{y) = 1, then 

Cov(r(y),X,/) = Cov(/(y),X,/) 

-Cov{g{F-iy)-PN^MF^{y),Vfix)),XJ) 

= Coy{fiy),XJ) - g{F-iy) - Pm^mF^ {y)^X^). 

The conclusion will therefore follow once we prove, for every y € M, 

Cov(/(y),X,./) = 5(F-(y)-P;v.MF"(y),^..), 

Coy{f{y)-p''{y)f{x),X.J) = g{F%y)-P^^MF%y),X^). 

As the two arguments are similar, we just prove the first equality. The 
map F^ can be decomposed as follows: 

F-{y) = PT^aM^F^iy) + PidM.F^iy), 

where 

i 

Pr^dM^F^y) = Y,Cov{f{y),X,f{x))Xj,,, 
i=i 

PT.aM.F''{y)= E Cov{f{y),X,f{x))Xj,,, 
j=i+i 

and the orthonormal basis {Xi^x, • ■ • , ^q,x) is chosen so that the set {Xi^x, ■ ■ ■ , ^i.x) 
forms an orthonormal basis for T^dMi and (Xi^i^x, ■ ■ ■ ,Xq,x) forms an or- 
thonormal basis for Tx dM-^, the orthogonal complement of Tx dMi in TxM. 

Further, because g{Xx,Vx) = for every Xx € TxdMi and Vx € NxM, it 
follows that 

PN.AiF^y) = PN^AiPiaM^F^y) 
and for every Xx G Tx dMi 

g{F^{y) - PN^MF'^{y),Xx)=g{PT^dM,F^iy),Xx) 

+ g{Pla^,f-{y)-PN^MPiaM,F-{y).Xx) 
= g{PT^dM,F^{y),Xx) 

i 

= J2Cov{fiy),X,f{x))g{Xj,x,Xx) 
i=i 

= Coy{f{y),Xxf). □ 

As noted above, the independence between V/|g^/-(x) and the process 
allows us to remove the conditioning on f\aMi{^) the expression for 
Diff whether / has constant variance or not. 
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Corollary 4.2. Suppose f is a Gaussian process satisfying the condi- 
tions of Theorem 3.3. Then, 



k 

Diff^,M(^)| < E / E(| det(-V2/|aA/„i?,(^))|lcf««] 



(32) i=o 



where 



Cf«« = {u< fix) - g{PT^9M.F^x),Vf{x)) < W{x)}. 
If f has constant variance, then F^[x) = S T^M and 

Cf^^ = {n</(x)<t^(x)}. 

Proof. The only thing that needs to be proven is that, in the event 
Cf^^ the condition {u < f{x) < W{x)} can be replaced with 

{u < fix) - g{PT^aM,F%x), V/(x)) < W{x)} 

and the conditioning can be removed. 

The reason that the above replacement is justified is that, on the set 

{V/|OM.(^)=0} 

f{x) = fix) - g{PT^9M.F^x),Vf{x)). 

Further, f{x)-g{PT^dMiF''{x),Vf{x)) is independent of VfiQM^{x), and 
Lemma 4.1 implies that W{x) is also independent of V f\dMiix)- Therefore, 
the conditioning on V fdhhix) can be removed. □ 

We are now ready to prove the following theorem. 

Theorem 4.3. Let f he a Gaussian process with constant, unit variance, 
on M and let f = f\M be such that f satisfies Assumption 2.1. Then, 

1 



liminf-n-^log|Diffo/(n)| > - 1 + 



where 



(33) <(/,x)= sup Var(r(y)) 

y€M\{x} 

and 

(34) a2(/)^supcT2(/,x). 
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Proof. We must find an upper bound for (32). Writing 

"^"^ f\dM,,Eii^) 

= V2/|aM„i?.(^) - E(V2/|9Af„S,(^)l/(x)) +E(V2/|9A^„iJ,(x)|/(x)) 

= vV|aA/.,i=;,(^) + /(^K-/(x)/ 

(cf. [22]), and applying Holder's inequality to (32) yields, for any conjugate 
exponents p, q, 

k 



Diff/,Af(n)| <^ / EE(/(x)'l{/W>n}) 



i=0 

xE(|detr,_,-(-V2/|aA/.,i5.(x)-/(x)/)n'/^ 

where detrfc(^) is the kth. "det-trace" of the square matrix A which is defined 
to be the sum of the determinants of all k x k principal minors of A. 
Define 

IJ, sup e( sup riy)). 

xeM \ y£M\x / 

For 

the Borell-Tsirelson inequality implies that 

P{W{x) >u)< 2e-("-^+)'/2'^'(/'^). 
Recalling that f^{x) = f{x), for such u, it also follows that 
E(/(^)'l{/(x)>.}) < C,^x^-ie-("-'^^)'/2. 
Putting these facts together, for any conjugate exponents p, q 
|Diff/,M(n)| 



k 

X " 



^ / j2^{\detv,^,{-W''f\9MMx)-f{x)I)ffPdn.{, 

• n J dMi ; r, 



The result now follows after noting that we can choose q close to 1, and 
u{q) so that, for u > u{q), the remaining terms are arbitrarily small loga- 
rithmically, compared to n^. □ 
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Theorem 4.3 provides a lower bound on the exponential decay of Diff j^a/(u). 
We believe that the lower bound is generally tight when a maximizer of cr^(/) 
occurs in dM^, in the sense that the term corresponding to dM^ in the sum 
defining DiS f^M{u) in (25) is exponentially of the same order as the upper 
bound; however, we were unable to prove this conjecture as it seems difficult 
to establish the sign of the error of the lower-order terms. In the piecewise 
smooth setting, it is therefore still open as to whether the liminf^^oo in 
Theorem 4.3 can be replaced with lim„_^oo as we cannot rule out the pos- 
sibility that some terms in the sum (25) cancel each other out, leading to 
a faster rate of exponential decay. Although we have not settled the issue 
completely, these situations seem somewhat pathological. 

5. Examples. In this section, we compute o"^(/) for some simple exam- 
ples, strengthening earlier results of [13, 16, 17]. Before turning to the exam- 
ples, however, we discuss the relation between cr'^{f) and the critical radius 
of a tube around M when / is assumed to be centered with unit variance. 
Specifically, we describe the geometry of the situation in the case of "global 
overlap," that is, when the supremum 

f^c(/)=sup sup Var(r(y)) 

x&M y£M\{x} 

is achieved at a pair {x*,y*), x* ^ y*. 

5.1. Geometric picture in the case of global overlap. Here, we describe 
the notion of "global overlap" and describe the geometry of the process / 
near pairs {x*,y*) achieving the critical variance a1{f). Roughly speaking, 
this situation occurs when M, the parameter space of /, "wraps around 
itself" and, for some x G M there is a point y G M that is close to x in 
the L^-metric but far in terms of geodesic distance from x. To describe the 
geometry involved in this situation we turn to spherical geometry in i/j, the 
RKHS of /. Recall that Hf is defined by the reproducing kernel condition 

{R{sr):R{t,-))jj^=R{t,s) 

and there exists an isometry that maps Hf onto the linear span Hf of 
{/(x),x € M} in L^(0,^, P). Without loss of generality, then, we can de- 
scribe the geometry in terms of Hf. Let ^' : M — Hj C L?'{il, J-, P) denote 
the map 

x ^ f{x). 

If S{Hf) is the unit sphere in Hf, and / is centered with unit variance, 
then ^(M) C S{Hf ), and our standing assumptions, namely that / is 
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and p{x,y) 7^ 1 if a; 7^ y, imply that ^ is a piecewise embedding. Further, 
it is not hard to see that 

so that the tangent space ry(2.)^(M) is spanned by {Xi^xf-,---iXk^xf) for 
some basis {Xi^^, ■ ■ ■ ,Xk,x} of T^M. 

We denote the orthogonal complement of Tj-(j.-)^'(Af) in T^(^x)Hf by r^^j^'(M), 
and the orthogonal complement of Tj(^.)^'(M) in Tf(^x)S{H f) by A^j(^)^(M). 

Given a point f{x) and a unit normal vector Vf(^x) ^ ^f(x)^{^) we denote 
the geodesic, in S{Hf), originating at f{x) in the direction Vf(^x) by '^f(x),Vf(^^) ■ 
That is, 

^/W.^'/w W = cost • /(x) + sint • < t < vr. 

As discussed in [21], up to a certain point along Cf(^x),Vff^^)i the points on 
cj(a;) metrically project uniquely to x. That is, for t small enough, the 
unique point on ^{M) closest to Cf(^x),vf(^){^) is f{x)- We denote the largest 
t for which this is true by 6{f{x),Vf(^x)) a-i^d we call it the local critical 
radius (angle) at f{x) in the direction Vf(^.j.)- Taking the supremum over all 
directions ^ '^i^f{x)^i^^)) obtain 6{f{x)), the local critical radius 
at f{x), 

(35) 0,{f{x))= inf 9{fix),Vf^x)) 
and the global critical radius 

(36) ec = e,{^{M))= inf 9,{f{x)). 

x£M 

The relation between these critical angles and a'^{f,x) in (33) is given in 
the following lemma. 

Lemma 5.1. Suppose f = f\M ^-^ the restriction of a centered, unit vari- 
ance Gaussian process on M , to M , a piecewise k-dimensional sub- 
manifold of M . Suppose that "^-.M—^Hf is a embedding. Then, for all 
X e dMk, 

al{f,x) = cot' {9,{x)). 
Proof. For x G dMk, N^M = T^M^ is a linear space, therefore 
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Furthermore, because / has constant variance F^{x) = 0. Putting these facts 
together shows that 

^ ^ f{y) - P^iy)f{x) -g{F-{y),Vf{x)) 
_ fjy) - P''{y)f{x) - Eii Cow{f{y),Xif{x))Xif{x) 

for some orthonormal frame field {Xi, . . . ,Xk) on M. 

Turning to the picture in terms of geodesies, fix f(x) and ^/(a;)- Suppose 
that for a certain t the point Cf(^x),Vf(^^) (^) does not metricahy project to f{x). 
This imphes there is a point f{y) € ^{M) such that 

dicf{x),Vf^,^it)Jiy))=cos~'^{{cf(^^)^^^.^^^{t),f{y))jj^,) 

< '^(c/(x),,;^(,) {t),f{x)) 

= COS"^((cj(^)^„^^^j it),fix))jj^). 

Alternatively, 

{cf{x),Vf(^^^it)Jiy))Hj = cost ■ {f{x), f{y))Hf + sint ■ {vf^^)J{y))jj^ 
= cost • p{x, y) + s\nt- {vf(^^),f{y))jj^ 
it), fix)) 

= cost. 

After a little rearranging, we see that this is true if and only if 

^^^^^(_!!fMif^y)hi^ 

Therefore, 

^/i / N N {vf{x),fiy))Hf 

cotec{f{x),Vf(^^))= sup — — 

yeM\{x} ^-Pix,y) 

Taking the supremum over all ^ ^f{x)^iM) we see that 

cot ^c(/(2;)) = sup sup ^ f{x),fiy))Hf 

^fM(^S{Nf^,^^{M)) yGM\{x} 1 - P{X, y) 
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where PNf^^^^{M) represents orthogonal projection onto A^j(2.)^'(M). There- 
fore, 

Hf 

cot e^{f{x))= sup 

y€M\{x} {^-P{x,y)r 

and it remains to show that 

\\PNf^^^^iM)f{y)\\Hf 

= Var (^f{y) - p{x, y)f{x) - CoY{f{y),XJ{x))XJ{x)^ . 

This, however, fohows from the fact that Nf(^,j.'j^(M) is the orthogonal com- 
plement (in Tf(^^-^Hf) of the subspace = span{/(x), . . . ,Xfc/(x)}, 
and the fact that 

k 

p{x, y)f{x) + Covifiy),Xifix))XJix) 
1=1 

k 

= {fix)Jiy))H,fix)+Y.{fiy),XJix))HfXJix) 

i=l 

= PLj{y). □ 

We now consider the case when M is a manifold without boundary and 
the supremum in 



(37) cot^(0c) = sup sup 



\PNf^^)-^{M)f{y)\?Hj 



x&My&M\{x} {^-p{x,y)Y 

is achieved at some {x*,y*). Thinking of the critical variance as the cotangent 
of some critical distance on M, then, geometrically, the tube of radius 6c 
should self-intersect along a geodesic from x* to y* in such a way that 
the tube viewed locally from the point x* shares a hyperplane with the 
tube viewed locally from the point y*. Alternatively, at the point of self- 
intersection the outward pointing unit normal vectors should be pointing in 
opposite directions. 

The simplest way of seeing this is to think of M as just two points {pi,P2} 
in M?. In this case the tube of radius r around M consists of the union of two 
discs of radius r. When r = d{pi,p2)/2 the two discs self-intersect at exactly 
one point, and the unit normal vectors are pointing in opposite directions. 

We can make this statement precise in the following lemma. 
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Lemma 5.2. Suppose {x*,y*) achieve the supremum in (37). Then, 
fiy*) = cos(20e) • fix*) + sm{2e,) ■ v}^^.^, 

where 

""^(^^^ " \\PN,,^.,nM)f{y*)\\H, 

is the direction of the geodesic between f{x*) and f{y*)- For any X^* € Tn^*M 
and any Xy* G Ty*M , 

Coy{X,, f,f{y*)) = {^,{X,*),f{y*))H,=0, 

Cov(X,./,/(x*)) = {^,{Xy,),f{x*))H,=0. 

Further, the partial map 

p^\y)^p{x\y) 
has a local maximum at y* , and the partial map 



py*{x) = p{x,y* 



has a local maximum at x* 



Proof. If {x*,y*) achieve the supremum in (37), then there exists a 
point z equidistant from f{x*) and f{y*) and unit vectors G Nj^^^t-^'ifiM) 

and Wf(^y*) € N f{^y*^'^{M) such that 

z = cos 6'c • f{x*) + sin 6'c • = cos 6'c • /(y*) + suiOcWf(^y*y 

It is not a priori obvious that w^f^y*^ G Nj(^y*-^^{M). However, if this were 
not the case, there exists t G M, t^y, such that d{t, z) < 9c- This contradicts 
the assumption that (x*,y*) achieve the supremum in (37) and the critical 
radius is 9c- The proof that this is a contradiction is left to the reader, 
though it follows similar lines to the argument at the end of this proof. 

To prove the first claim, therefore, it is enough to show that the unit 
tangent vectors c/(^*),^^j^.j (6'c) and c/(y*),^„^(^.J (6'c) at z satisfy 

u = c/(x*),«;^(,.) (Oc) + c/{r ),«;^(,.) i^c) = G T,Hf- 
Suppose then that u^O- A simple calculation shows that 

cos{29c)-{fix*),f{y*))H, 



{u,f{x*))Hf = {uJ{y*))Hf 



sm( 



If u 7^ 0, then f{x*),f{y*) and z are not on the same geodesic and the 
triangle inequality implies that 



cos(20c)<(/(:r*),/(y*))^ 
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which imphes 

{uJ{x*))H, = {u,f{y*))H,<0. 
Consider the geodesic originating at z in the direction u* = u/\\u\\f{f 
Cz,u* (t) = cos t • z + sin t ■ u*. 

Assuming u* 7^ 0, 

{c,,u*{0),f{y*))H, = {cz,u'{0)J{x*))hj = {u*J{x*))<0. 
This imphes that for sufficiently smah \s\, s < 0, 

{f{x*),Cz,u'-{s))Hf = {f{y*),Cz,u*{s))Hf >COs6c. 

For such an s, there exist distinct points x{s) and y{s) such that 

dicz,u* is),x{s)) < 9c, d{cz^u* is),y{s)) < 0c 

such that the geodesic connecting x{s) and Cz^u*{s) is normal to M at x(s), 
and the geodesic connecting y{s) and c^^„*(s) is normal to M at y{s). 
Without loss of generality, we assume that 

d{cz,u-'{s),x{s)) < d{cz,u*is),y{s)). 

In this case, the geodesic connecting y to Cz^u*{s) is no longer a minimizer 
of distance once it passes the point Cz^u*{s), which is of distance strictly less 
than 6c from y{s). That is, for some point z(s) along the geodesic connecting 
y{s) and Cz^u*{s), beyond Cz^u*{s) but of distance strictly less than 9c from 
y(s), there exist points in M strictly closer to z(s) than y{s). We therefore 
have a contradiction, as the critical radius of ^(M) is 9c. 

To prove the second claim, we note that f{y*) is a linear combination 
of f{x*) and '^'/(x*) which are perpendicular to every vector in rj(j,.)\I'(M). 
Similarly, f{x*) is a linear combination of f{y*) and Wf[y*) which are per- 
pendicular to every vector in Tj(^.*)^'(M). The fact that the partial maps 
p^* (y) and p^* (x) have local maxima follows from the same contradiction 
argument used above. □ 



5.2. Centered stationary processes on [0,T]. In this section f{x),x G 
[0,T], is assumed to be a centered stationary process with unit vari- 
ance and covariance function R. As in [13] we change the time scale so that 
Var(/(x)) = -ij(0) = l. 

Fix a point x £ (0,T), in which case the process is given by 

p. ^ ^ fiy)-Rix-y)f{x)-R{x-y)f{x) 
^ l-R{x-y) 
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and the critical variance at x is given by 

2,,^^_ _ Var(/(|/)|/(x),/(x)) 



c^c(/>a;)= sup 

j;g[o,T]\{x} {^-R{x-y)) 



l-R(tf-R{tf 

~ -.T^<T-. (1 - R{t)Y ■ 

The critical variance in the interior is 

(38) a2(/,(0,r))= sup al{f,x)= sup ^ p7,^fi^^' • 

xG(0,T) 0<t<T {^-J^yt)) 

A local version of critical variance (j^ ^^^(ir, (0, X*)) is obtained by letting 
t^O in (38), that is, 

^g.loc(^,(0,T)) = lhn ^ =R^'\^) - 1 

where the last conclusion follows from some simple calculus. An alternative 
interpretation of the quantity [qj.(x, (0, r)), is the following: 

a,%,(x,(0,r))=Var(/(x)|/(rE)). 

The local critical variance at the end points {0,T} needs slightly more 
attention. We consider the point x = without loss of generality. The normal 
cone at X = is 

^o[0,T]=Uc;^croM. 

c<0 

For ease of notation, we just write A'^o for Ao[0,T]. 
The projection onto Aq is just 

As 

Cov(/(0),/(y)) = -i?(y) 
the process f^{y) is given by 

_ /(y) - i?(2/)/(0) - l{^(^)<o^i?(y)/(0) 



r\y) 



{l-R{y)? 

The local critical variance at x = is given by 

^ l-R{ty-R{t)^ + max{R{t),Oy 

O-c(j>0)= sup 

0<t<T {i-R[t)y 

(39) 

1 - Rjt)^ - R{tf 
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The inequality above implies that (T^(/, 0) > (T^(/, y), < y < T. Therefore 
the critical variance cr^(/) is attained at the end points t = 0,T. 

Note that this does not exclude the case that the critical radius is also 
attained in the interior < t < T. Under some circumstances, the critical 
variance is attained everywhere, as demonstrated in the following lemma. 

Lemma 5.3. Suppose f is a centered, unit variance stationary pro- 
cess on M such that —R{0) = 1 and R{t) < for all t>0. Then 

^c(/) = <ioc(/,^) = Var(/(x)|/(x)). 

Proof. As noted above, we just have to compute cj^(/, 0). Because 
R{t) < we see that 

2..^^_ l-Ritf-Ritf 



C7,(/,0)= sup . 

o<t<T [i^-R{t)) 

Now suppose that this supremum is achieved. Lemma 5.2 implies that for 
any t that achieves this supremum, R{t) = 0, and t is a local maximum of 
R{t). Let 

T* = {t> 0\R{t) = 0, R{t) = cos(26'c)}. 
As T* is a closed set there exists a minimum value of T*: 

r = minT* >0. 

Because t* is a local maximum of R, there exists e > such that 
R{t)<R{t*) yt£{t* -e,t*). 
But R{t) is assumed to be nonincreasing and we have 

R{t)>R{t*) yt<t*. 

Therefore 

R{t) = R{t*) yt£{t* -e,t*). 
This, however, contradicts the minimality of t*. □ 



Remark. Presumably, the same ideas as in [5] could be used to gen- 
eralize Lemma 5.3 to a larger class of covariance functions, though we will 
leave that as an exercise for the reader. 
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5.3. Isotropic fields on with monotone covariance functions. In this 
section, we compute the critical variance for centered, unit variance isotropic 
processes restricted to a compact, convex set M. In particular, we prove 
the following proposition, which shows that the exponential behavior of 
Diff f^M{u) for such processes is determined solely by the conditional variance 
of the second derivative, given the field. 

Proposition 5.4. Let f be an isotropic process on M™' that induces 
the standard metric on satisfying Assumption 2.1 and let f = f\j^ he the 

restriction of f to a compact, convex set M with piecewise smooth boundary. 
If the covariance function 

R{\\x\\) = CovU{x + t)J{t)) 

is monotone nonincreasing, then the critical variance is attained locally and 
is given by 

^2 . 



^c(/)=Var(^(0)|/(0)). 



Proof. Fix t G M° . Similar computations to those in Section 5.2 show 
that 

Because R is assumed monotone nonincreasing, the arguments used in Lemma 5.3 
imply that 

sup Var(/*(s))=Varf^(0)/(0)' 

We therefore turn to the boundary dAI, assumed to be piecewise smooth. 

Fix t G dM and s^t. Let Xf be the unit vector in TiW^ in the direction 
s — t. Since M is convex, € 5jM and 

Xt \\F\s)eStM, 

that is, Xf is parallel to F*, so that Pn^mF^{s) =0. Therefore, 

We see that we do not incur any additional penalty for the local critical 
radius at the boundary points if M is convex. Again, the arguments of 
Lemma 5.3 imply that 

^2. 



a2(/,t) = Var(^|^(0)|/(0)) 



The conclusion now follows. □ 
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